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HALL EQUILIBRIA: SOLUTIONS WITH TOROIDAL AND POLOIDAL 
MAGNETIC FIELDS IN NEUTRON STAR CRUSTS 

K.N. Gourgouliatos, 1,2 A. Cumming, 1 M. Lyutikov 3 ' 4 and A. Reisenegger 5 

RESUMEN 
ABSTRACT 

We present Hall equilibrium solutions for neutron stars crusts containing toroidal and poloidal magnetic field. 
Some simple cases are solved analytically while more complicated configurations are found numerically through 
a Gauss-Seidel elliptic partial differential equation solver. 
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1. NEUTRON STAR CRUST 

The crust is the upper layer of a neutron star 
(~lkm) and consists of a highly conducting crystal 
lattice. Lorentz forces exerted by the magnetic field 
are balanced by the elastic crust and the field evolves 
through Hall effect ( Goldreich &; Reisenegger|[l992 ). 



2. HALL EVOLUTION 

In the Hall effect, the magnetic field is advected 
by the electric current, which is carried by electrons. 
Ohmic diffusion, for neutron star crust applications, 
is much slower and is treated as a secondary effect. 
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In a system dominated by Hall drift, equilibrium 
occurs when the Hall term is equal to zero. The 
field might settle into this state after several Hall 
timescales. We solve for Hall equilibrium inside the 
neutron star crust while requiring that the field con- 
nects to an external vacuum dipole field. 
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3. ANALYTICAL SOLUTIONS 

We write the magnetic field in terms of the 
poloidal flux function W and the poloidal current 
cI/2: 



B = V* x V</> + /V<; 



(2) 



Neglecting the Ohmic term and demanding that 
<9tB = 0, we obtain the Grad-Shafranov equation. 
Defining the Grad-Shafranov operator A* = d 2 + 
Sln 2 e dg (-^a)) Hall equilibria are given by 



A** + II' + r 2 sin 2 6n e S' = . 



(3) 



where a prime denotes differentiation with respect 
to \P. / and S are arbitrary functions of Analyt- 
ical solutions which connect smoothly to a vacuum 
field outside the star are possible for linear choices 
of S oc \& and 1 = 0. These fields have the angular 
structure of a dipole and correspond to uniform ro- 



tation of the electron fluid (Cumming et al. 20041 



Analytical solutions for fields with mixed poloidal 
and toroidal components fully confined within the 
star are possible for / cx ^ and bear similarities with 



previously known MHD solutions (Gourgouliatos et 
aL pHOl ), see figure [lj 



4. NUMERICAL SOLUTIONS 

To avoid the drastic simplifications of analytical 
solutions, we numerically solve the Grad-Shafranov 
equation. The azimuthal field is confined to closed 
tori as the external vacuum cannot support currents. 
The numerical scheme reproduces the analytical re- 
sults and has some similarities with the poloidal 



fields found in Hall simulations (Pons & Geppert 



2007 Kojima & Kisaka 
plotted in figures |2|and|3) 



2012). The solutions are 
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Fig. 1. Meridional sections of the star showing surfaces 
of constant W which are also the poloidal field lines. Left: 
The purely poloidal field solution corresponding to S oc 
W and 7 = 0. Right: The fully confined mixed poloidal 
and toroidal field corresponding S oc W and / oc W. In 
both plots n e =const. and the thickness of the crust is 
0.2, assuming a stellar radius of unity. 
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Fig. 2. Plot of the magnetic flux function at the equator 
W(r, it/2) as a function of the radius of the star found 
numerically. / is chosen to be / = a(W — Wo) 1 ' 1 where 
Wo is the value at the equator, S' =const., the thickness 
of the crust is 0.2 and n e = 1 — r 2 . The solid line is 
the analytical solution, the black crosses is the numerical 
solution for the same case (a = 0) while the red circles 
and the blue stars correspond to a = 50 and a = 100 
respectively. 




Fig. 3. Plots of the poloidal field lines with toroidal field. 
The colour represents the value of W. The left plot corre- 
sponds to a — 15 and the right to a = 45. The toroidal 
field is hosted in the enclosed red loops. The choices of 
parameters are as in figure [2] 

5. NEUTRON STAR MAGNETIC EVOLUTION 

The solutions presented, subject to their stabil- 
ity, may represent long-term states in the magnetic 
evolution of neutron stars. Despite their similari- 
ties, they do not coincide with MHD equilibria, thus 
a phase transition during the formation of a neutron 
star from a fluid to a solid state will be followed 
by Hall evolution that may lead to Hall equilib- 
rium. The Hall equilibrium will be slowly decaying 
by Ohmic diffusion from the start. Eventually, when 
this effect has already caused a substantial reduction 
of the field strength, the Ohmic term will dominate 
over the Hall term, and then the field structure will 
be essentially that of the lowest-order Ohmic mode 
rather than the Hall equilibrium. A more extended 
and in-depth discussion of this work will be pub- 
lished soon (Gourgouliatos et al., in preparation). 
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